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ABSTRACT 


Hydrodynamic loads induced on large underwater structures 
by impulsive motion may be significant design factors. Such 
loads may be induced by earthquake excitation or may result 
from acceleration produced while lifting an object in the 
sea. 

The theoretical approach to the calculation of these 
loads is outlined and numerical results are presented for 
several submerged configurations including a vertical cone, 
a sphere, and a vertical circular cylinder. 

Numerical results for these submerged structures are 
presented in the form of a dimensionless load parameter or 
йесі mass coefficient. Results corresponding to a number 
of different water depths are presented to show the rather 
ле effect of the relative water depth on the hydro- 
dynamic force. It is shown that the effect of the free 
water surface is to reduce the hydrodynamic loads in 
comparison to the corresponding infinite depth values. 

Experimental results obtained by vibration testing are 
presented for a submerged sphere, a cone and a vertical 
cylinder. These results show excellent agreement with the 


theoretical results. 
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NOMENCLATURE 

Description 
characteristic dimension of the body 
added mass coefficient 
depth of submergence 
source strength function 
Шене соло 
acceleration of gravity 
Green's Function 
depth of water 


modified Bessel function of the second kind 
of order nm ` 


distance of line of action above bottom 
moment vectom 

unit normal vector 

dynamic pressure, dimensionless 

dynamic pressure 

fluüid velocity vector 

surface of body 

time 


velocity vector of object 


= coordinates, dimensionless 


c^&/g 


= coordinates of a point on the surface of 


the object, dimensionless 


fluid density 





symbol Description 


С = frequency 

$ = velocity potential, dlmenslonless 
$ - velocity potential 

5 = angular vector velocity of body 
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I. INTRODUCTION 


When an object as accelerated Through at Ве 
are, in general, two types of forces that are recognized, 
one being a drag component and a second, the inertial 
component. The nature of the flow produced by the time 
dependent motion of an immersed rigid object and the relative 
contribution of these two components of force is generally 
Переа to be strongly dependent on the amplitude of the 
relative fluid motion in comparison to the characteristic 
lineal dimension of the object. For example, Keulegan and 
Carpenter [Ref. 1] found that, for harmonic motion of the 
fluid past a. fixed circular cylinder, major flow separation 
did not occur and the forces were well represented by 
potential flow values provided the amplitude of the motion 
was less than about a half diameter. If the amplitude 
is small enough the fluid does not move in one direction 
far enough for separation and wake development to occur 
and, accordingly, an inviscid flow analysis is a valid 
mathematical model. | 

Accelerated motion Of a solid body ша а fluid occurs 
Rar example when a bottom mounted structure is excited by 
an earthquake. Other applications include the accelera- 
tion of ships and accelerated motions of objects being 
lifted in the ocean environment. Southerland [Ref. 2] has 
described the effect of added mass on the mechanics of 


lifting such loads as a submarine rescue chamber. 
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John [Ref. 3] has shown that ene velocity potential 
associated with tino Wem e motion of a rigid DeC 
immersed in a fluid with a free surface may be represented 
by a surface distribution of point wave sources for arbi- 
trary values of frequency. The source distribution (or 
weighting) function is obtained in this formulation from 
the solution = integral equation resulting from the 
kinematic boundary condition applied on the immersed surface. 
Using this technique, numerical results have been obtained, 
for example, for a semiellipsoid oscillating in a free sur- 
face in water of infinite depth by Kim [Ref. 4], the 
intended application being to surface ships. Numerical 
results for a bottom mounted hemisphere have been presented 
by Garrison and Seetharama Rao [Ref. 5] and results for 
several axisymmetric configurations by Milgram and Halkyard 
[Ref. 6]. Garrison and Chow [Ref. 7] have recently presented 
both theoretical and experimental results for wave forces 
acting on a fixed submerged oil storage tank. 

In the work cited, the hydrodynamic force coefficients 
are represented as a function of a dimensionless frequency 
parameter (or equivalently, a dimensionless wave length 
parameter) of the form v = oa where o denotes the 
frequency of the motion, a s the characteristic lineal 
dimension of the submerged object and g denotes the gravi- 


tational constant. This dimensionless parameter is essentially 


a form of the Froude number. In all cases, however, the 


TIT 





theory has a fundamental limitation. The source potential 
oscillates with a wave length inversely proportional to 
the dimensionless frequency of the motion v, and consequently, 
a reasonable partition T of the immersed surface limits 
the validity of the numerical scheme to small to moderate 
Paes of frequency. Specifically, 1 has been found 
using the Haskind's relations [Ref. 8] and an energy balance 
for purposes of checking accuracy, that results can be 
obtained with reasonable accuracy up to about v = 2.0. 
Beyond this, the accuracy of the numerical results becomes 
questionable. 

Although the range v = 0 to about 2.0 includes most 
applications in ship hydrodynamics and wave forces acting 
on large structures, the Wmequencies encountered, for 
example, in earthquake excitation corresponds to rather 
large values of v in the case of large structures. Accord- 
ingly, the development of an alternate method which includes 
both the effect of the free surface and bottom on the 
hydrodynamic forces and is valid for large values of the 
dimensionless frequency v is the primary objective of the 
present thesis. Therefore an attempt has been made herein 
to obtain numerical results for several different objects 
for the asymptotic case of very large values of v, i.e., 
Por the =antinite Hroude number case, 

High frequency added mass — for two-dimensional 
shapes oscillating on a free surface have been presented 


in a series of papers by Landweber et al. [Ref. 9, 10, 1, 
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and 12]. The motivation for Emus work was to provide 
hydrodynamic coefficients for two-dimensional ship forms. 

For three-dimensional objects the only results known 
are those of Waugh and Ellis [Ref. 13] for the case of a 
sphere in infinite depth water and Jacobsen [Ref. 14] 
more tne case or the cylindrical pier in various water 
depths. In the former case, the theoretical method consisted 
of using successive image doublets to account for the effect 
of the free surface. However, the case of vertical motion 
only was studied and their method was limited to a sphere. 
In the latter case, the solution was represented by a series 
of Bessel functions and is applicable to a vertical circular 
cylinder, extending between the bottom and free surface. 
However, e numerical computations were made using tables 
of Bessel functions with Price arguments which have 
proven misleading. In fact, Jacobsen's numerical results 
are in error in several instances and therefore, his 
analysis as well as the corrected formulae are presented 
herein. 

In this thesis the theoretical formulation of the 
impulsive hydrodynamic forces acting on a rigid object ої 
arbitrary shape is presented. The effect of the bottom 
and free surface is taken into consideration and the force 
coefficient is shown to depend on a frequency of oscillation 
parameter. The problem is first formulated for arbitrary 
values of the frequency and the asymptotic form of the 
solution for both large and small values of the frequer .y 


is discussed. A computer method for numerical evaluation 


13 


of the force coefficient for submerged objects of arbitrary 
shape is discussed T" numerical results are presented 

for a number of specific geometric shapes including a 
circular cylinder, cone, and sphere. High frequency 
experimental results are presented for a vertical circular 
мәнде, Sphere, and cone. These experimental results 


are compared with the theoretical results. 
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II. THEORETICAL CONSIDERATIONS 


A. BOUNDARY VALUE PROBLEM 

The geometry involved in the problem under consideration 
is shown in Figure 1. A rigid object of arbitrary shape 
having characteristic lineal dimension a is submerged to a 
depth d in water of depth h and caused to move in a small 
amplitude harmonic motion. The instantaneous motion of the 
rigid object is specified by the lineal velocity Ü and 
angular velocity Q. The object may intersect the bottom 
ОТ free surface. 

Bene an unseparated, incompressible and irrotational 


—- 


flow, a velocity potential (x,t) defined by 


=> 


а = Vo (1) 
=> 
where q denotes the fluid velocity vector, may be introduced 


and must satisfy the continuity equation, 


v^8 = 0 (2) 


within the fluid region. The bars over the spatial variables 
> 


denote dimensional quantities and x = ime + JEA + KX. 
On the rigid bottom surface described by X5 nen 
fluid velocity in the vertical direction must vanish and, 


accordingly, the corresponding kinematic boundary condition 


on the velocity potential is 


39 (x, »—h, xX, N) = 
9X» 


0 (3) 


ЕР 


The velocity of a point on the immersed surface is given 


by 


r3| + 


(H) 


> 
where P denotes the position vector as shown in Fig. l. The 


appropriate kinematic and dynamic boundary condition on the 
NNIEPSedosumbace. therefore, takes the form 
> 


art) vn (5) 


an 


where n = Gun = iia + Kn, denotes the unit normal vector on 


3 
mae surface of the object and is directed outward inte the 
Fil jd region. 
Assuming the amplitude of the motion of the object to 
be small, the kinematic and dynamic boundary condition applied 


at the free surface may be linearized to give the well-known 


free surface boundary condition: 


E — D = 
39 (хуз0,Х3»б) З 13 9(x,,0,x5,t) m (6) 
9х2 9t 


Furthermore, if the object oscillates with frequency o, the 


time dependence of the velocity potential may be separated 


-ict 


and written as ReL6 (x4 ,x5,x.)e ] where Re denotes the real 


part. In which case Eq. (6) becomes 


2— 
ОЕ Ona _ | 
3x 1t 573' -O 7x $(x4,0,x5) 7 0 (7) 





where the unbarred quantities represent the coordinates ` ade 


dimensionless with a. 
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Fig. 1 - Definition of Geometry 


In general it is apparent from Eq. (7) that the solution 
to the na value problem established must depend on the 
dimensionless frequency parameter о ИЕ. However, there are 
two extreme cases where the results become independent of the 
frequency. These correspond to the case where c*a/g > œ jn 


which case Eq. (7) becomes 
$ (x4 ,0,x-) = 0 · (8) 


For the other extreme where ave + 0 the free surface appears 


as a rigid boundary as Eq. (7) becomes 


39: (x1,0,x2) = O (9) 
9х 
2 . 

If consideration is restricted to either of the boundary 
conditions Eq. (8) or Eq. (9) the velocity potential may be 


written as the sum 


=> 


MR ve) = W, (t) $4 Q0, i a UD (10) 


where each term in the sum represents the potential associated 
Deich the ee component of body motion. (i = 1,2,3 refers to 
linear motion components in the X1»X55X4 directions, 
respectively, and i = 4,5,6 refers to the angular motion 


components about the X »X55X5 axes, respectively.) The time 
dependent function W is defined as 
W,(t) = U,(t) a 
i = 42.3 Chi) 
= —2 
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in which the symbols U; and Q. denote the components of the 


i 
lineal and angular velocity vectors of the body, U and on 
respectively. The potentials $4 are dimensionless quantities. 
It is important to note that the application of either of 
the two boundary conditions, Eqs. (8) and (9), makes the 
problem homogeneous in time. Consequently, the assumption 
of harmonic motion may be relaxed and the functions W, (t) may 
be allowed to be arbitrary functions. 
Кере с етан На (10) into Eqs. (2), (3), (5) emcee or 


(9) and introducing the dimensionless variables defined as 


x q ae р-н d= d/a (12) 


the boundary value problem for фі (1 - 1,2,...0) пауйбе 


specified as follows: 


v^ $4 04,4, x4) = 0 


26. 
$5(x1,0,x,) 7 0, (o9) or A (4:0) 


9X, 


364 Gt 7h, x3) 


x 0 (13 a-d) 
ОСО) 
— 172° 3 _ m (on the immersed surface) 
an i 
where m, =n, , i = 1,2,3 
ту = (xətd)n. - ха 
"p | Ио 


пб = хуп, - (x, 4 d)n, 
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(14 a-d) 





The boundary conditions (13b) represent two different 
possibilities on the gue surface, the first being the high 
frequency condition and the second the rigid boundary (or zero 
frequency) condition. In as much as the method of solution 
involving either boundary condition is similar, both solutions 


are developed simultaneously in the following. 


B. GREEN'S FUNCTION SOLUTION 

The solution to the boundary value problem (13) may be 
carried out by use of a Green's function. The Green's 
function represents the potential for a point source and 
these sources are distributed over the immersed surface 
ccording to the distribution function ey: According to 


this concept the potential $ may be written as 


00 = & Jf г, (2) atas (15) 


where (ГЕ?) denotes a point on the immersed surface, G denotes 
the Green's function and dS = a8/a* denotes the surface area 
element made dimensionless with the characteristic body 
dimension a. 

In order that Eq. (15) represent a solution to Eqs. (13) 


it is necessary that G satisfy the equation 


v^Q(X;2) = 6(%-?) (16) 


where 6(x) denotes the Dirac delta function, as well as 


boundary conditions (13c) and (13b). 
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The Green's function satisfying Eq. (16) as well as 
Eq. (13b) and (13c) is obtained by use of successive images 


of sources as depicted in Fig. 2 and is given by 


00%) -2- R £25 261)" m +=) + C (= + yl 


n ТІ n n 


(1) 


The corresponding form of G satisfying the rigid boundary 
condition Eq. (13e) (o = 0 оп X, = 0) is also obtained by 
successive images and is similar to Eq. (17) except that all 
Ae eSources have positive strengths. For this базе, che 


Green's function takes the form 


б(Х;Е) = š + = + > | = + = + = + =) (18) 
1 n=l 2 3 Y 
n n n n 


where in both Eq. (17) and (18) 


ш = V(x, - D + (x, - з)" 
В = ios t (X, = LAS 


2 
1 0 + (x, + Е. 


R. = V we + (x, - 2nh - p (19) 
pc MANT Ex 
Э 2 2 


ел) 
II 


ye 


za 


4h+n(+) 


2h-n(*) 





Note: 

When the free surface acts 

as a rigid boundary all 

imaqe SOUrCes are Ositive, 
-(2h-n) (-) 4 J P 


. mus 412% S 4- 
те. в r Lage source ОООО 





Ry = \ „© + (x, + 2nh - 


n 


Nee) 
5 2 2 
n 
Since G satisfies Eq. (13a-d) the potential ф, as given by 
Eq. (15) must also satisfy these conditions. Thus, it remains 
to select the function f such that Eq. (13d) is satisfied. 
Application of this boundary condition yields the following 


Fredholm integral equation of the second kind 
> ii > CE) 5 
X; de Е 
-fi(x) t Em JJ г, (2) SA з dS = 2m, (x) ; зов о 5M UE 


which is to be solved for e This integral equation is to 
be satisfied at all points X on the immersed surface. 

Because of the form of G, Eq. (20) 1s rather complex and 
it is therefore necessary to carry out the solution to the 
integral equation numerically. Proceeding in se direction 
the immersed surface may be partitioned into N area elements 


of size 4S, and the integral in Eq. (20) written as the sum 


=f. +f, %,, = 2m >, ЖІ 52 (21) 


where i,j = 1,2,...N and n refers to the six possible modes 

of motion of the object. The transformation from the integral 
equation to a summation is possible since f is a well-benened 
Manctlon. It follows from Eq. (20) that the coefficient matrix 
is defined by 
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й af f задо) 
а) "о 15, m i as (22) 


The function 3G/an in Eq. (22) is obtained in a straight- 
forward manner by differentiation of Eq. (17) or (18). 


Once the coefficient matrix a is obtained from Eq. (22) 


ij 
a standard digital computer subroutine may be used to invert 
the matrix equation, Eq. (21), to obtain the source strength 
f at points on the immersed surface. For purposes of evalua- 
tion of G and 3G/3n by use of (17) and (18) it was found that 


about 15 terms were adequate. The data presented herein, 


however, were obtained using 30 terms. 


C. HYDRODYNAMIC FORCES AND MOMENTS 
The dynamic fluid pressure is obtained from the linearized 


form of Bernoulli's equation as 


Q2 
e 


(25) 


FJ 
! 
! 
т 
Q2 
ct 


where the velocity squared term has been neglected. 


Corresponding to the ¿e component of body motion, the pressure 


is obtained by use of Eq. (10) and (23) as 


P (no sum) AO (24) 


The pressure may further be written in dimensionless coeffi- 


cient form as 


D) Ж 
=-9,(x) , 12 1,2,...6 (25) 





Ee 
p, (x) p . 
oW 


F 
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The hydrodynamic force acting on the object is obtained 
by carrying out the integration of the pressure over the 
immersed surface. Using Eq. (24) the force vector associated 


with the gen component of body motion is 
B. = =? ||. P.ndS, er 5 (26) 
с 9 


The corresponding result for the moment vector is 


Й, = 3 fT Р. їх В 8, ј = 1,2,...6 (27) 
9 а J 


Using Eq. (24) as well as the definition of m. as given in 


Eq. (14), the qub component of force associated with the 


ЈЕ" component of acceleration is given in dimensionless 


coefficient form as 


Aig x s m; dos 
тыл 


and the corresponding expression for moment is 


ow = bm; dS, i=} (29) 
J pa № i j Tee ane 6 


where w is defined by Eq. (11). 


DS (28) 
2 б 


Cae 
|| 


| 
~ 
л 
` 
СА 


Equations (28) and (29) define the inertia (or added mass) 
tensor having, in general, 36 elements. The dimensionless 
Poefficients Ass defined by Eq. (28) represent the dimension- 
less force coefficient and Eq. (29) represents the dimension- 


less moment coefficient. These coefficients are generally 


referred to as added mass (or moment of inertia) coeff. ents. 
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Although these coefficients. are made dimensionless by use of 
a in the definitions (28) and (29), there are cases, depending 
on the body geometry, where other length scales are more 
appropriate. It is also standard practice in many cases to 
use the displaced fluid volume. However, in all numerical 
results presented herein the definition of Ај is specified 
on the figure so that although the definitions vary, no 
Ке пелоп should result. 

For purposes of numerical evaluation of the coefficients 
given in Eqs. (28) and (29) the integrals may be written as 
the summation, where the integral sign has been replaced by 


the sum on the k index. 


Ay; = (64m). AS), › (30) 


However, in order to apply Eq. (30) it is necessary to first 
evaluate Ф; at all of the nodal points on the immersed 
surface. For this purpose the surface integration in Eq. (15) 


is also converted to the summation 


Gy m ACD, о K (31) 


where the coefficient matrix is defined as 


ва = JUN G(X,;2) as | (32) 
9, 5 


Once the coefficient matrix is evaluated by Eq. (32), 


the summation indicated in Eq. (31) may be carried out t 
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determine $; at all nodal points on the immersed surface. 


The coefficients A,, are then obtained by use of Eq. (30). 


J 
Although there are 36 possible added mass coefficients 

defined by Eqs. (28) and (29), many are zero on account of 

Symmetry and most off-diagonal elements are small. Also, 


even for a body of arbitrary shape As; is symmetrical so that 


J 
at most only 21 of the coefficients need be evaluated. That 


d 
is, Eq. (13e) along with Eq. (28) and (29) gives 


9%. 
As y e $m. dS "E 93 va (33) 


By use of the Green's theorem applied to the fluid region 


external to the body, ít can be shown that 


TM . 39; 
КО жи э Ж Оше > (31) 


since on the free surface either ¢ = 0 or 9ф/39п = 0 (depending 
on the boundary condition being enforced) and 3¢/dn =.0 on 
the bottom. Therefore, in view of Eq. (34) and (33) it is 


apparent that the added mass tensor is symmetrical, 


A (35) 


pee JACOBSEN'S SOLUTION 
At this juncture it is appropriate to comment on the 
analytical solution presented by Jacobsen [Ref. 2] for the 


Special case of a vertical circular cylinder extending between 


at 





a rigid bottom and a free surface Ри SOlution pro еле 
the only known analytical result for the problem under 
consideration and, therefore, is of considerable value as a 
point of comparison for the numerical results presented 
herein. However, it was found that a number of errors were 
present in Jacobsen's work and it was necessary to re-develop 
Ea re-evaluate the equations. The expression for the added 
Mass coefficient associated with horizontal acceleration in 
perms of the notation of the — paper and in terms of 
Bessel functions having real as opposed to imaginary argu- 


ments appears as 


ТОМЕ | 1 nT 
hy = 3 6 | DES K, (Sh) (36) 
П=1 59.5 





where A41 denotes the added mass coefficient made dimension- 


tees With the displaced volume, the factor В, is аеГїпеа ас 


? . 
_ ТУП 
хер + к бр > 


TA I 


and Ко апа K; denote modified Bessel functions of the second 
kind of order zero and one, respectively. This expressicn 
agrees with Jacobsen's corresponding expression with the 
exception of a factor of 2.0 and the type of Bessel functions 
involved. It was found that more serious sign errors occurred 


in Jacobsen's expression for the line of action of the hori- 


zontal hydrodynamic force and, accordingly, his numerical 


28 


results are also in error. . The correct expression for the 
elevation of the line of action (or center of added mass) in 


terms of the notation of the present paper is 


l nti 2 nT 
| IC pl BK, (sp) 


n=1,3 





% = (38) 


пт 


5 
1 
n BAK (sp) 


n=1,53,9 


where 2 = 1/Һ, & being the elevation of the line of action 
above the bottom and h the water depth. 

For completeness the expression for the dimensionless 
иге оп the surface of the cylinder is given in terms of 


the notation of the present paper as 


> ВЕ 

сов 6 (-1) “ соз| ШТ (nty)] B K. (Z) (39) 
E gy 1 1 

T n=1,3,5 

Euere 8 denotes the angle measured with respect to the x-axis 
in the x-z plane. 

Equations (36) and (37) were programmed for numerical 
evaluation on a digital computer so that an adequate number 
of terms could be evaluated without concern over numerical 
error. These results are plotted in Figs. 8 and 9 for 


comparison with the source distribution results developed 


23 


herein. It was noted that numerical results obtained from 


figures presented by Jacobsen [Ref. 2] are in error when 


compared to numerical results obtained from Eqs. (36) and 


(38). 
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ТТТ. EXPERIMENTAL APPROACH 


Experiments were conducted to determine the added mass 
coefficient of a circular cylinder, a cone and a sphere in 
proximity of a free surface and rigid bottom. The basic 
experimental approach consisted of mounting the test 
object on flexible beams and measuring its natural frequency 
Sou vibration in air and water. The reduction in natural 
frequency when placed in water was then related to the 
added mass of the water. 

The test rig was used to determine the added mass 
coefficient corresponding to the horizontal motion is shown 
in Fig. 3 and 4. The deis was mounted in this test rig 
on a thin 6 inch long by 2 inch chord vertical strut. Two 
large steel channels were placed perpendicular to each 
other and welded to a 1/2-inch steel plate to form a rigid 
support for the spring-mass system. These channels were 
Supported at four points on the top edge of a shallow, water- 
Echt box (24" x 70" x 48"). The cylinder, cone, and sphere, 
mounted on a pair of flexible beams, were bolted to this 
supporting structure as shown in Fig. 3. (The cylinder 
is shown.) 

Four strain gages were mounted on the upper end of 
one of the flexible beams and connected into a Wheatstone 
bridge. The bridge was then connected to a carrier pre- 


amplifier and the frequency of the bridge output was measured 
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Fig. 3 - Experimental equipment. 
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by use of a Monsanto programmable counter-timer and in some 
cases by the use of a Visicorder. A ten second duration 
was set on the counter-timer and, depending on the depth of 
water and test object, the readings ranged from about 110 
to 180 cycles. The time duration programmed on the counter 
was highly accurate so that the major error in the process 
mas considered to result from reading the number of cycles 
occurring in the 10 second interval to integer values. 
However, since the starting of the interval was random, 
Baas error was minimized by repeating each test run six 
times and recording the average reading. It was noted that 
mnie spread in readings from the counter was never greater 
than one cycle in the ten second interval for a given 
@oniiguration. 

A calibration curve for the system was developed by 
plucking the mass (test object) in the air with the 
addition of pre-measured weights attached to the test 
object. As a result it was unnecessary to rely on the 
theoretical natural frequency equation for the spring- 
mass system. In general, the test procedure was considered 
to be highly accurate and good repeatability was obtained 
in all cases. The calibration curves are shown in Appendix 
E 

The same procedure was used for the determination of the 
added mass coefficients for the sphere in vertical motion 
but a different apparatus and tank was employed. The 


parallel beam mount which carried the rod with attached 
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inch diameter plexiglass sphere is shown in Fig. 5 and 
6. The tank used in this experiment was 4.0 ft. deep and 
represented, for practical purposes, an infinite depth 


case. The sphere was oscillated in the vertical direction. 
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Fig. 5 - Experimental apparatus used in heave tests. 
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IV. DISCUSSION OF RESULTS 


A. CIRCULAR CYLINDER 

Theoretical and experimental values of added mass 
coefficients have been generated for a number of configurations 
and are plotted in Figs. 7-20 and tabulated in Tables l- 
18. These include results for a vertical circular cylinder, 
a sphere, a hemispherical bottom-mounted vessel, and a 
ea object. 

The theoretical and experimental results for the added 
mass coefficient for a vertical circular cylinder are shown 
in Figs. 7-12 and Tables 1-7. The calibration curves are 
presented in Fig. T. Two different cylinders of the same 
overall dimensions and materials were used. They produced 
calibration curves which were slightly different -resuits 
Since their natural frequencies of vibration were different 
аш ретт? positions in the test rie could not осе "Ор али 
exactly. 

The added mass coefficient corresponding to the case 
of the horizontal acceleration as based on the cylinder 
radius cubed is presented in Fig. 8, both with a free sur- 
face and with a rigid upper boundary. Results for two 
different grid sizes are shown, 120 and 192 effective nodal 
Points, and iittle difference was found to exist. “fl сап 
be seen that the Beiter results for the "rigid boundary" 
case closely follow the straight line relationship (th) 


representing the classical two-dimensional result. The 
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results corresponding to the "free surface" case show that 
the оғ the free surface Is Lo redueor pora deze 
coefficient and for larger values of h the results appear 
сере closely approximated by the straight line relation- 


Ship A (rh - 3.0). 


je. 

The same type of result in a slightiy different form is 
presented in Fig. Y. In this figure the added mass coeffi- 
cient made dimensionless in the more common manner, by use 
of the displace volume ra“h , is plotted for two different 
grid sizes along with Jacobsen's results, Eq. (36). Also, 
the experimental results obtained from vibration testing 
are presented for comparison. 

The results show that the effect of the free surface 
is to reduce the added mass coefficients and agreement 
between the experimental and theoretical results is.:shown 
ГО De excellent. On the same figure the computer results 
for the rigid boundary case are presented. The classical 
result for two-dimensional flow gives А] = 1.0 and, 
therefore, the deviation from this value gives the percentage 
Afmerteal error in ¿he computer results directly. Тһе 
figure indicates a maximum of about 3.0 per cent error zt 
the e Tow values of h. 

There is no doubt thart 22 Of the computer results 
could be improved by means other than the most obvious: 
increasing the number of area eem ES However, any gain 


in precision is generally offset by increased computer 


storage requirements, complexity and computer run time. The 
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present method is considered to be quite simple and versa- 
tile ноте Е to configuration and yet possesses adequate 
accuracy for most purposes. 

For Une case sor horizontal motion, the sine oi action 
of the horizontal force is also a parameter of interest. 
When the upper boundary acts as a "wall", the force acts 
at the center, i.e., at 2 = .5. However, for the free 
сасе савзе Еле pressures are reduced near the free зоти 
amd the line of action is somewhat below center. huis re nu 
is presented in Fig. 10 as obtained using 252 nodal points 
and is compared with Jacobsen's corrected result, Eq. (38). 
The figure SNOWS Che agreement Lo ве excellent. 

A second degree of freedom is also possible for the 
vertical cylinder. The cylinder may be allowed to rotate 
about some convenient axis with angular acceleration. In 
which case an added moment of inertia may be defined. 
Numerical results have been generated, using 192 nodal points, 
for this coefficient for rotation about an axis passing 
through the center of the cylinder at its base. These 
theoretical results are presented in Fig. 11. The results 
show that the free surface has the effect of reducing the 
hydrodynamic reaction when compared to the rigid boundary 
case. This is the expected result since the free surface 
tends to reduce the pressure differential near the top ОП 
the cylinder. | 

Finallg, in order to show how the maximum pressure on 


the surface of the cylinder varies with depth, Fig. le as 
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been prepared. In this figure the dimensionless pressure 
coefficient as defined by Eq. (25) is plotted for the case 
of horizontal acceleration. It may be noted that the 
classical solution corresponding to the rigid boundary case 
yields the maximum pressure of p, = 1.0 over the complete 
length. Pig. 11 shows that the pressure corresponding to 
the free surface case begins at zero and appears to approach 
Шейоу а5 ће distance from @he free surface and depth ineredse, 
This is the expected trend since the effect of the free 
surface diminishes at large distances. Eq. (39) is plotted 
on the same figure for comparison and the agreement appears 


cor be quite good. 


БІ SPHERE 

Results for a sphere accelerating both horizontally and 
vertically are shown in Figs. 13-17 and tables 8-14. The 
calibration data is presented in Fig. 13 (vertical) and 
Fig. 14 (horizontal). For the case of vertical acceleration 
in deep water, Fig. 15 shows experimental results obtained 
by vibration testing for comparison with theory and the 
agreement appears to be quite good. It may be noted that 
the effect of the free surface is always to decrease the 
added mass coefficients while the proximity of a rigid 
boundary has the opposite effect. All of the computer results 
for the sphere have been generated using 264 nodal points. 

The well-known classical solution for a sphere gives an 


added mass coefficient of .5 for the case of a fluid of 
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infinite extent. Accordingly, as d becomes large, the 
Coefficient plotted in Pig. 12 3s expected ров рр о на 

this closed form result but approximately a 4.0 per cent 
номтеь-оп 125 noted. This error is due to numerical Ina ours 
acies and approximations used in the computer program. 

The same kind of results are plotted in Fig. 16 as 
generated by the computer, the only difference being the 
finite depth, h = 4.0. The same general trends as in Fig. 
12 are observed when the sphere is near the free surface; 
the rigid boundary results are greater than .5 while the 
free surface results are less than .5. However, as the depth 
is increased and the sphere approaches the bottom, the 
added mass coefficient increases in either case. 

Theoretical and experimental results for the added 
mass coefficient .in enis ons acceleration corresponding 
to a bottom mounted sphere is shown in Fig. 17. The 
experimental results were obtained by use of the test rig 
shown in Figs. 3 and !. The sphere was mounted on the lower 
plate by use of a 6.0 inch strut and the sphere was adjusted 
рО within .02 inches of the tank floor. It may be noted 
that the results show good agreement with theoretical 
results and have the same general trends as previous results 


except for the break off in A for the rigid boundary case 


11 
as h becomes small. This reversal in trend is’ caused by 
the unwetting of the sphere occurring at values of h less 


than 2.0. 
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С. HEMISPHERE 

Figure 18 and Table 15 present the added mass coeffi- 
cient associated with a bottom mounted hemispherical object 
as a function of water depth for both the rigid boundary and 
free surface case. The results appear somewhat similar 
to previous ones with the rigid boundary having the effect 
of increasing the added mass coefficient and the free surface 


having the opposite effect. 


D. CONE 

Conically shaped bottom mounted, surface piercing 
structures have been proposed for use in the Arctic on 
account of the ice-breaking capability of the shape. This 
structure shown in Fig. 19 is a conically shaped tower 
similar to those which have been proposed as drilling rig 
ои тептабјопг for deployment in the Arctic. The Computer 
results were generated using 192 effective nodal points 
distributed over the immersed surface. The tabular data 
is listed is Tables 16-18. The calibration curve is pre- 
sented as Fig. 20. | 

The theoretical and experimental added mass coefficient 
(based on actual displaced volume) for the conically shaped 
tower is presented in Fig. 19 as a function of water depth 
and these results are in good agreement. The results for 
the rigid boundary case shows coefficients which are much 
гелет than those corresponding со сре tree SUriace OC E 


Moreover, for the free surface case the line of action f 
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the hydrodynamic force was found to actually lie below the 
bottom. This rather unexpected result is caused by the 

fact that the pressures on the lower portion contribute to a 
moment opposite that induced by pressures on the upper part. 
In general, as'the object is accelerated horizontally, 
positive pressures are produced on the front side while 
negative pressures occur on the back side. Therefore, since 
the pressure as well as area which contributes to a negative 
moment is much greater than those tending to produce a 
positive moment, the result is a moment tending to tip 

the top of the cone in the direction of the acceleration. 
This moment is equivalent to applying the horizontal force 


at a distance t below the bottom. 
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V. CONCLUSIONS 


A practical numerical analysis has been developed using 
three dimensional sources and image sources to represent 
the flow produced during impulsive motion of rigid struc: 
‘Cures of arbitrary shape immersed in water in proximity to 
a free surface. The bottom and free surface have been taken 
into account and it is shown that the proximity of the free 
surface has a sizable effect and always tends to reduce the 
added mass coefficients in comparison to the rigid boundary 
ease. The experimental results obtained by vibration 
pe Scing for awspnhere, Cylinder, and a cone were found to 
agree well with the theoretical results and the results for 
the eylinder are in agreement with analytical results. 
The results are directly applicable to such practical 
problems as the determination of earthquake loading of rigid 
submerged and semi-submerged bottom mounted structures as 
well as the calculation of required forces to lift objects 


of arbitrary shape. 
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APPENDIX A 


TABLES OF COMPUTER GENERATED DATA 


TABLE 1 


CYLINDER: COMPUTER GENERATED DATA (Horizontal Motion) 
192 Effective Points 


A /о0, nach 


"e zog 


Free Surface 
Boundary Cond. 


0.910 
0.896 
0.879 
0.861 
0.057. 
0.807 
0.761 
0.699 
85225 
0.512 
0.402 
0.254 


A JoU nach 


Jam H 


oodd 
Boundary Cond. 


1.000 
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-001 
.004 
ED 
. 014 


ШЕ 


=. e He P HKH ы 


3027 
1051 
12027 
1.062 
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h Ату 
Free Surface 
Boundary Cond. 





q = 


0.472 
0.469 
0.465 
0.461 
0.455 
0.449 
0.422 
0.434 
0.425 
0.419 
0.414 
0.412 
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TABLE 2 


CYLINDER: COMPUTER GENERATED DATA (Horizontal Motion) 
120 Effective Points 


mo іі + ал с ~ CO КО 
сео о 


оны m 


мл 


о м 


A, y=F,,/0U,12R А 1=2 1/20. та? в = 2 

Prec асе Solid Free onoi en 

Boundary Cond. Boundary Cond. Boundary Cond. 

0.935 1.006 0.476 

0.908 1.002 0.473 
0.889 1.000 0.469 
0.869 1.000 0.465 
0.844 120903 0.459 
0.813 12007 0.453 
0.769 1.014 0.446 
0.704 1.922 0.437 
0.597 т 0.427 
0.515 1.036 0.421 
0.404 1.045 0.416 
02255 1.069 0.410 


5 





TABLE 3 


CYLINDER: COMPUTER GENERATED DATA (Horizontal Motion) 
252 Effective Points 
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0.907 0.999 0.471 
0.893 0.998 0.468 
0.877 0.998 0.464 
0.858 0.999 0.460 
0.835 1.002 0.455 
0.804 1.006 0.449 
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TABLE 6 


CYLINDER: COMPUTER GENERATED DATA (Rotational Motion 
about vertical axis) 192 Effective Points 


Acg7Mec/ (02. 18^h^) Acg7Mgg/ (08,78*h?) 
Free Surface Solid 
Boundary Cond. Boundary Cond. 

0.266 0.322 

0.250 0:320 

0.248 0.319 

0.237 0.317 

0.225 0.314 

0.209 ооо 

0.190 0.308 

0.166 0.302 

E 0.294 

0.110 0.288 

0.083 0.282 
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TABLE 7 


CYLINDER: COMPUTER GENERATED DATA (Pressure 
Coefficients, Free Surface Boundary 





Condition) 
ED 
РІ ^ 23018 а = 2/8 h = ћу 
0.929 ДЕО 2 5.0 
0.925 4.167 5.0 
0.916 3.611 5.0 
0.900 3.056 5.0 
0.875 2.500 5.0 
0.835 1.944 5.0 
0.769 1.389 o 
0.649 0.833 5.0 
0.394 орао 5.0 
2.833 Олон 3.0 
2.500 0.841 3.0 
РОТ 0.827 3.0 
1.833 0.804 3.0 
1.500 0.769 5.0 
1.167 OS 3.0 
0.833 0.637 30 
0.500 2-50 Se 
Gener, 0.279 зоб 
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TABLE 12 


COMPUTER GENERATED DATA (Horizontal Motion, 


Bottom Mounted ) 


A, TF. (pU, 1582/3) 


Free Surface 
Boundary Cond. 


0.669 
0.667 
07565 
0.648 
07625 
0.555 
07529 
0.440 
Ose 
0.292 
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TABLE 15 


HEMISPHERE: COMPUTER GENERATED DATA (Horizontal Motion) 


E Á, ,7F44/ (00, 2182/3 A, 744 / (00, 2182/3) 
Free Snri ace Soliq 
Boundary Cond. Boundary Cond. 

6.0 0.516 0.518 

3.08 0.515 0.519. 

4.0 0.514 0.520 

DO 0.510 0.526 

а 95 0.505 Олова 

1.7 | 0.481 0.567 

1.4 0.455 : 0.610 

Шо 0.421 0.674 

О 0.366 0.901 
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- TABLE 16 


CONE: COMPUTER GENERATED DATA (Horizontal Motion) 


B АР 1/0019 А11=211/2 0% 
Басе асе Solid 
Boundary Cond. Boundary Cond. 
on | 0.885 0.901 
4.0 0.869 0.922 
0 0.845 0.960 
2.0 0.766 1.109 
1.5 0.664 1.405 
К 25 0.578 1. 7 
h | $ Ё 
Free Surface Solid 
Boundary Cond. Boundary Cond. 
6.0 0.336 0.336 
4.0 0.338 0.335 
3.0 0.55) 0.555 
2,0 0.348 0.324 
1-5 0.370 0.305 
1,25 0,394 0.281 
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CYLINDER: 


Frequency (HZ) 


18. 
Ша- 
15. 
15. 
1. 
14. 
По. 
ie. 
ШЕ! 
12. 
Lae 
11. 


№0. 
19. 
15. 
14. 
14. 
15: 
ie 
mec. 
T2: 
I: 


00 
15 
70 
30 
DE 
15 
85 
80 
25 
70 
80 
10 


25 
63 
27 
86 
10 
80 
48 
47 
04 
67 


APPENDIX B 


EXPERIMENTAL DATA 


TABLE 4 


CALIBRATION DATA (Horizontal Motion) 


Cylinder (A) 


Cylinder (B) 


. 67 


Added 


о Ооо > P P O O 


TT EU nn мю F° > © o 


Weight (lbs) 


250 

. 793 
Оса 
ea 
. 814 
.043 
.296 
‚508 
. [90 
.329 
. 360 
.388 


0 

. 861 
SEI 
. 866 
. 351 
a 
.880 
.888 
. 373 
. 890 


о Oo O F FP гг NN WWW fF LF VI 


NU E eS MI 


.30 
. 80 
„20 
ТЕ 
15 
ето 
„25 
.00 
„65 
„25 
„90 
YO - 
.00 


250 
.80 
720 
O 
15 
.75 


CYLINDER: 


TABLE 5 


EXPERIMENTAL DATA (Horizontal Motion) 


Cylinder 


Frequency (HZ) 


12: 
12. 
15. 
33% 
14. 


14 


(A) 


10 
50 
00 
48 
10 


о 
15. 
192 
16. 
162 
Pt 
17. 
or 


30 
66 
10 
[5 
20 
78 
00 


Cylinder 


12 


. 43 
12. 
13. 
15: 
14. 
15. 


80 
40 
2 
60 
OF 
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(B) 


A 


O O O ооо о о о о о со с O 


CO. C CO CON tO 


Added Wt. 
11 Displaced Wt. 


. 835 
oo 
. 804 
115 
чо 
"оса 
.618 
1505 
‚548 
„441 
. 390 
221 
‚000 


SOUS 
шб 
MS 
„135 
. 694 
.679 





TABLE 9 
SPHERE: CALIBRATION DATA (Vertical Motion) 


Frequency (HZ) Added Weight (1lbs.) 


19. 
19. 
18. 
10% 
17. 
16. 
16. 
16. 
15. 
15. 
15. 
14. 
14. 


N N HM . FP ым о о о о 
O — хл to) O 21 хл м» O NN UL 9 о 


no 


= Оо о ы = о = -і пло Оо о о 
го 
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TABLE 10 


SPHERE: EXPERIMENTAL RESULTS mx Motion, Infinite fluid 
Depth 


h д _ = Added Weight 
22 Displaced Weight 


6.0 0.514 
Б.О 0.480 
4.0 0.514 
Sn 0.504 
20 0.500 
LS 0.475 
о 0.463 
1525 0.416 
1.0 0.30 
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TABLE 13 


SPHERE: CALIBRATION DATA (Horizontal Motion, Bottom Mounted) 


Frequency (HZ) e Added Weight (1bs.) 
14.65 0.379 
14.35 0.592 
14.00 0.828 
13.95 0.864 
ПЕ 65 1.077 
13.40 1.313 
15.00 1.625 
22 [5 7.0753 
12.25 | | 2.390 
11.85 | 2.875 
11.45 3. 390 
20.75 4.405 
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TABLE 14 


SPHERE: EXPERIMENTAL DATA (Horizontal Motion, Bottom Mounted) 


. added wt. 
h L; ^11 = displaced wt. 


3.212 12.05 0.644 
3.40 12.00 0.661 
EIS 12.05 0.644 
2.90 Too 0.624 
2.67 12.15 0.612 
2.50 12.20 0.605 
2.33 12.25 0.555 
2.17 lo s 0.563 
2:00 12.45 0.536 
1.83 O 0.465 
o o 12090 0.421 
1.50 lett 0.355 
1.33 13.55 0.287 
I ЕНЕР 0.230 
1.00 14,15 0.179 
0.00 15.32 0.00 


T2 


TABLE 17 


CONE: CALIBRATION DATA (Horizontal Motion) 


Frequency (Hz) Added Weight (lbs.) 
17.61 0.418 
16.89 0.631 
16.50 0.867 
15053 1. 332 
PST 1.665 
14.62 1.912 


ШЕ 





© © © © O 


TABLE 18 


CONE: EXPERIMENTAL DATA (Horizontal Motion) 


Frequency (Hz) Ај) = ЕТ 
Wt. 
MU E 0.281 
.64 16.03 0.214 
48 16.89 0.142 
33 17.79 0.095 
.11 18.66 0.065 
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